#include "substitution_matrix.h"

double substitution_matrix_distance (unsigned char a, unsigned char b) {
  double ret;
  if(a < 24 && b < 24) {
    ret = MAX_BLOSUM_SCORE-blosum62[a][b];
  } else if( a == b) {
    ret = 0.0;
  } else {
    ret = 20.0;
  }
  //  printf("%d, %d = %f\n", a, b, ret);
  return ret;
}

